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BACKGROUND OF THE INVENTION 

1 Field of the Invention 

This invention relates to computer-assisted methods for identifying and designing small 
molecule mimetics of erythropoietin. 

2 Description of Related Art 

Erythropoietin (EPO) is the primary regulator of the proliferation and differentiation of 
immature erythroid cells. EPO is produced in the fetal liver and in the adult kidney » 
response to hypoxia (low oxygen levels in blood or tissue). It circulates in the blood 
stream where it targets the EPO receptor (EPOR) on committed progenitor cells in the 
bone marrow and other hematopoietic tissues. Recombinant human erythropoietin 
(rHuEPO) is widely used in therapy of patients with anaemia due to chronic renal failure, 
cancer chemotherapy and AZT treatment 

^m^M^ ^om^^^^^^^ indudes receptOTS 

for other hematopoietic growth factors such as interleuldns (ILs). colony stimulating 
factors (CSFs) as well as growth hormone prolactin and ciliary neurotrophic factor 
(CNTF). The structural architecture of this family of receptors consists of three modules: 
a Ugand binding extracellular domain, a short trans membrane region and a large 
cytoplasmic domain. It has been proposed that the extracellular domain of this 
supcrfanfly comprises two discrete domains each containing approximately 100 residues 
that fold into a sandwich consisting of 7 antiparallel P-strands with the topology of an Ig 
constant domain. Members of the family share two characteristic motifs in their 
extracellular domain: a pair of conserved disulfide bridges in the N-terminal domain, and 
a WSXWS box (where X is any amino acid residue) in the C-terminal domain. For most 
members of this receptor superfamily. oUgomerization of one or more polypeptide chains 
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is essential for forming high affinity receptor complexes. A homodimer complex has been 
demonstrated to be the active fonn of hGHR and a similar model has been suggested for 
G-CSF, prolactin and EPO receptors. 

Erythropoietin induces dimerization of two EPO receptor molecules, which results in 
subsequent phosphorylation of the cytoplasmic domains by the association with two 
tyrosine kinase (JAK2) molecules to initiate a cascade of events that leads to the relevant 

biological. 

Given the importance of erythropoietin, it would be very desirable to be able to identify 
molecules capable of binding the EPO receptor and eliciting the response normally 
elicited by EPO. 
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SUMMARY OF THE INVENTION 

The invention features methods for identifying molecules which will bind to the EPO 
receptor and act as a EPO mimetic. Preferred EPO mimetics identified using the method 
of the invention act as agonists of the EPO receptor in one or more in vitro or in vivo 
biological assays of EPO activity. Preferred mimetics are molecules lacking peptide 
bonds, i.c, are non peptidic mimetics. Preferred peptide mimetics have 15 or fewer, more 
preferably 10 or fewer amino acids. 

The methods of the invention entail identification and design of molecules having a 
particular structure. The methods rely on the use of precise structural information derived 
from x-ray cry^allographic studies of the extracellular domain of EPO receptor (amino 
acids 1 to 225) complexed with a peptide, EMP1 (EPO Mimetic Peptide 1; described 
below), which acts as an EPO mimetic. This crystallographic data permits the identifica- 
tion of atoms in the peptide mimetic that are important for EPO receptor binding and 
dimerization. More importantly, this data defines a three dimensional array of the 
important contact atoms. Other molecules which include a portion in which the atoms 
have a similar three dimensional arrangement similar to some or all of these contact 
atoms are likely to be capable of acting as an EPO mimetic. Moreover, one can use the 
structural information to design or identify molecules having even more EPO activity 
than the peptide mimetic described herein. 

The details of the preferred embodiment of the present invention are set forth in the 
accompanying drawings and the description below. Once the details of the invention are 
known, numerous additional innovations and changes will become obvious to one skilled 
in the art. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a flowchart showing a first method for identifying potential mimetics of 
erythropoietin using a computer system. 

FIGURE 2 is a flowchart showing a second method for identifying potential mimetics of 
erythropoietin using a computer system. 

Like reference numbers and designations in the various drawings indicate like elements. 
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DETAILED DESCRIPTION OF THE INVENTION 

Throughout this description, the preferred embodiment and examples shown should be 
considered as exemplars, rather than as limitations on the present invents. 

Described below is the crystal structure of a small peptide mimetic of EPO bound to an 
extraceUular portion of the EPO receptor. The peptide, EMP1 
(GGTYSCHFGPLTWVCKPQGG; SEQ ID NO: 1). is characterized by an mtramolecular 
Several lines of evidence suggest that EMP1 can act as an EPO mimetic. 
ToTeWe- EMPl competes with EPO in receptor binding assays and induces ceDular 
pr oliferation of cell lines engineered to be responsive to EPO. Both EPO and pep.de 
Luce a similar cascade of phosphorylation events and cell cycle progress^ m EPO 
responsive cell, Furhter, EMP1 demonstrates significant erythropoietic effects m rruce 
as Litored by two different in ^ assays of nascent red Wood cell production . Tta 
data, when combined, strongly supports the notion that the peptide hgand, wluch has a 
sequence unrelated to that of EPO. is capable of binding to and inducing an agorust 
conformation or assembly of EPO receptor. 

Design of small molecule miraetics 

TnestructureoftheEMPl dimer demonstrates that a molecule substantially smaller than 
natural hormone can act as an agonist and induce the appropriate biological response. 
The peptide is assumed to have a substantially smaller contact interface with the receptor 
than its natural hormone. The binding detenninanU m the EPO rector form an almost 
flat surface which is mainly hydrophobic in nature, without any cavities or charged 
residues that may help in design of a small molecule ligand to interact with the receptor. 

T« s simplified framework of interactions revealed by the structural data presented herein 
can be used to identify additional EPO mimetic, The atoms of EMP1 which are 
iro portant for binding to the EPO receptor and forming dimeric EPO receptor include 
thlinvorvedinthe contact between the EMP1 (peptide) and EBP (EPO receptor) and 
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those involved in contacts between the two EMP1 molecules in the dimeric complex 
(peptide-peptide contacts). In addition to the contacts listed in Table 2, the following 
EMP1-EMP1 hydrophobic contacts are significant: Tyr p \ Qy?"*, Phe M , Trp PlJ ,and Gyi" 15 
in each peptide. The following EMP |-EBP hydrophobic interactions are also'significant: 
Tyr w , Phe M , and Trp p,J in each peptide. It will be understood by those skilled in the art 
that not all of the atoms present in a significant contact residue need be present in a 
mimetic. In fact, it is only those few atoms which actually from important contacts with 
the EPO receptor which are likely to be important for mimetic activity. Those skilled in 
the art will be able to identify these important atoms based on the model of the dimeric 
EMP1-EPO complex which can be constructed using the structural data herein. 



Preferred mimetics will include atoms at postions similar to those of the EPO receptor 
contacting atoms of EMP1. Even more preferred mimetics will be structurally similar to 
the dimer of EMP1 found in the structure described below. This is because the 
dimerization of EMP 1 is an important factor in the diemerization of the EPO receptor. 

The methods of the invention employ a computer-based methods for identifying 
compounds having a desired structure. More JS p.pcffieaUy i ^ 
dimensional coordinates— of a„ subset of the atoms" in -the —peptide 
GGTYSCHFGPLTWVCKPQGG when the peptide is co-crystallized with a portion of 
the erythropoietin receptor comprising amino acids 1 to 225 of the human EPO receptor, 
to determine peptide and non-peptide mimetic candidates by means of computer methods. 

These computer-based methods fen into two broad classes: database methods and de novo 
design methods. In database methods the compound of interest is compared to all 
compounds present in a database of chemical structures and compounds whose structure 
is in some way similar to the compound of interest are identified The structures in the 
database are based on either experimental data, generated by NMR or x-ray crystallogra- 
phy, or modeled three-dimensional structures based on two-dimensional {i.e.. sequence) 
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data. In de novo design methods, models of compounds whose structure is in some way 
similar to the compound of interest are generated by a computer program using 
information derived from known structures, e.g. t data generated by x-ray crystallography 
and/or theoretical rules. Such design methods can build a compound having a desired 
structure in either an atom-by-atom manner or by assembling stored small molecular 
fragments. 

The success of both database and de novo methods in identifying compounds with 
activities similar to the compound of interest depends on the identification of the 
functionally relevant portion of the compound of interest. For drugs, the functionally 
relevant portion is referred to a pharmacophore. A pharmacophore then is an arrangement 
of structural features and functional groups important for biological activity, e.g., EPO 
activity. 

Not all identified compounds having the desired pharmacophore will act as an EPO 
mimetic. The actual activity can be finally determined only by measuring the activity of 
the compound in relevant biological assays. However, the methods of the invention are 
extremely valuable because they can be used to greatly reduce the number of compounds 
which must be tested to identify an actual mimetic. 

Dimerization of the EPO receptor is important for activity. Accordingly, preferred 
mimetics will be based on the structure of the EMP1 dimer as it is bound to the EPO 
receptor dimer. Thus, preferred mimetics have include important contacts from both of 
the RWJ 61233 peptides present in the structure described below. Such mimetics will 
favor dimerization of the EPO receptor. 

Programs suitable for generating predicted three-dimensional structures from two- 
dimensional data include: Concord (Tripos Associated, St. Louis, MO), 3-D Builder 
(Chemical Design Ltd., Oxford, U.K.), Catalyst (Bio-CAD Corp., Mountain View, CA), 
and Daylight (Abbott Laboratories, Abbott Park, IL). 
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*«*grams^it^te^^ 
bearingra desiredphanna^ 

S>i>xl/3DB"Unity (TriposAssodat^Sri^uisTKiO)^ 

Programs suitable for pharmacophore selection and design include: DISCO (Abbott 
Laboratories, Abbott Park, EL), Catalyst (Bio-CAD Corp., Mountain View, CA), and 
ChemDBSOD (Chemical Design Ltd., Oxford, U.K.). 

Databases of chemical structures are available from Cambridge Crystallographic Data 
Centre (Cambridge, U.K.) and Chemical Abstracts Service (Columbus, OH). 

De novo design programs include Ludi (Biosym Technologies Inc., San Diego, CA) and 
Aladdin (Daylight Chemical Information Systems, Irvine CA). 

Those skilled in the art will recognizefthat^theid^gn^ require slight 

structural alteration or adjustment of a chemical structure designed or identified using the 
methods of the invention. 

In general, chemical compounds identified or designed using the methods of 'the invention 
can be sythesized chemically and then tested for EPO activity using any of the methods 
described below. The methods of the invention are particularly useful because they can 
be used to greatly decrease the number potential mimetics which must be screened for 
EPO activity. 

The invention may be implemented in hardware or software, or a combination of both. 
However, preferably, the invention is implemented in Computer; programs executing on 
programmable computers each comprising a processor, a data storage system (including 
volatile and non-volatile memory and/or storage elements), at least one input device, and 
at least one output device. Program code is applied to input data to perform the functions 
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described above and generate output information. The output information is applied to 
one or more output devices, in known fashion. The computer may be. for example, a 
personal computer, microcomputer, or workstation of conventional design. 

Each program is preferably implemented in a high level procedural or object oriented 
5 programming language to communicate with a computer system. However, the programs 
can be implemented in assembly or machine language, if desired. In any case, the 
language may be a compiled or interpreted language. 

Each such computer program is preferably stored on a storage media or device (e.g., 
ROM or magnetic diskette) readable by a general or special purpose programmable 

,o ' computer, fcrconfi^ 

is read by the computer to perform the procedures described herein. The inventive system 
may also be considered to be implemented as a computer-readable storage medium, 
configured with a computer program, where the storage medium so configured causes a 
computer to operate in a specific and predefined manner to perform the functions 

is described herein. 

FIGURE 1 is a flowchart showing a first method for identifying potential mimetics of 
erythropoietin using a computer system. The method uses a programmed computer 
comprising a processor, a data storage system, at least one input device, and at least one 
output device, and comprises the steps of: 
20 (1) inputting into the programmed computer through an input device data comprising 
the three-dimensional coordinates of a subset of the atoms in the peptide 
GGTYSCHFGPLTWVCKPQGO when the peptide is co-crystallized with a 
portion of the erythropoietin receptor comprising amino acids 1 to 225 of the 
receptor, thereby generating a criteria data set (STEP 100); 
as (2) comparing, using the processor, the criteria data set to a computer database of 
chemical structures stored in the computer data storage system (STEP 102); 
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(3) selecting from the database; using a program suitable for searching three- 
dimensional databases to identify molecules bearing a desired pharmacophore 
(such as those described above or equivalents), chemical structures having a 
portion that is structurally similar to the criteria data set (STEP 104); 

(4) outputting to an output device the selected chemical structures having a portion 
similar to the criteria data set (STEP 106). 

FIGURE 2 is a flowchart showing a second method for identifying potential mimetics of 
erythropoietin using a computer system. The method uses a programmed computer 
comprising a processor, a data storage system, at least one input device, and at least one 
output device, and comprises the steps of: 

(1) inputting into the programmed computer through an input device data comprising 
the three-dimensional coordinates of a subset of the atoms in the peptide 
GGTYSCHFGPLTWVCKPQGG when the peptide is co-crystallized with a 
portion of the erythropoietin receptor comprising amino acids 1 to 225 of the 
receptor, thereby generating a criteria data set (STEP 200); 

(2) constructing, using a program suitable for generating chemical structure models 
(such as those described above or equivalents), a model of a chemical structure 
having a portion that is structurally similar to the criteria data set (STEP 202); 

(3) outputting to the output device the constructed model (STEP 204). 

Confirmation of Biological Activity 

In order to determine whether a molecule identified using the methods of the 
invention can act as an EPO mimetic, one or more in vitro or in vivo assays of EPO 
activity should be performed. For example, mimetic molecules should be able to 
stimulate proliferation of TF-1 cells (Kkamura et al. f J. Cell Physiol. 140:323, 1985) or 
B6Sut cells (Greenberger et aL, Proc. Natl. Acad. Sci. USA 80:2931, 1983), but 
preferably do not stimulate proliferation of cells which do not bear the EPO receptor. 
Thus, preferred rnimetics do not stimulate proliferation of Mo7e cells (Avanzi et aL, Br. 
J. Haematol. 69:359, 1988). 
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Potential nineties can also be tested in a murine model of erythropoiesu. In tbs assay 
a potential mimetic is administered to normal mice which express endogenous basal 
,evels of EPO. Reticulocytes are counted, preferably by flow cytometry, to determine 
whether the candidate mimetic increases reticulocyte levels. An increase in reticulocyte 
,evels indicates that the candidate mimetic is stimulating erythropoiesis. Because the 
riceused in this assay already express EPO, this assay may be relatively insensitive. As 
an alternative, candidate mimetics can be assayed in the -hypoxic-polycythemic mouse 
bioassay In this assay polycythemia is induced by conditioning mice in a hypobanc 
chamber to reduce endogenous EPO levels. Apotential EPO mimetic can be adminis- 
tered to a conditioned mouse. Incorporation of »Fe into blood serves as a measure of 
erythropoiesis. This erythropoiesis can be attributed to the candi^te m,rneuc. 

The assays described above are examples of soluble assay, Other assays for EPO 
activity known to those skilled in the art are also useful. 

In order to determine the biological activity of a candidate mimetic it is preferable to 
rae asure biological activity at several concentrations of candidate mimetic. The activity 
at a given concentration of candidate mimetic can be compared to the acuvity of EPO 
itself. 

Structural Data . ^ ... 

The coordinates for amino acids 1 to 225 of the human EPO receptor bound to peptide 
EMP, are presented in the attached appendix in standardarooKnayen database format 
Aisomdudedinthisappen^ These ordinate, 

can be used in the design and identification of EPO mimetics according to the methods 
of the invention. 

Structure of EBP-EMP1 Complex 

The extracellular fragment of human EPO receptor (EPO bmdmg protem, EBP) 
consfcting of residues 1-225. was expressed in Zstenti^ and purified as desenbed 
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(Johnson ct aL f Protein Express. Purif. 7:104, 1996). Rhomboidal-shaped crystals of an 
EBP complex with EMP1 were obtained in orthorhombic space group P2 I 2,2 I> with cell 
parameters a=59.2A, b=75.5A, c=132.2A, with two EBP and two peptide molecules in 
the asymmetric unit and a V H =2.8 AVdalton (Matthews, J. Mol. Biol. 33:491, 1968). The 
crystal structure was determined by multiple isoraorphous replacement (MIR) using two 
heavy atom derivatives (Table 1). Residues 1-2 and 19-20 of each peptide as well as 
residues 1-9, 21-23, 164-166, 221-225 of receptor molecule I, and residues 1-9, 21-23, 
133-135, 221-225 of receptor molecule II had poor or no electron density and are 
excluded from the structure analyses 

An important break in the electron density that affects the structure interpretation occurs 
for the three residues (Arg» -Gly» -Pro*) that link the amino terminal a-heiix to the first 
P-strand in Dl of both receptor molecules. A molecular packing diagram shows the 
proximity of a second non-crystallographically related dimer in the crystal that gives two 
possibilities of how this three-residue linker may be connected. The current choice of 
linker connectivity is based on a structure of another independent EBP-peptide complex 
at higher resolution (2.5 A), which shares a similar molecular packing, but for which the 
electron density is clear for these three residues. At present there are no experimental 
data to verify whether this N-terminal cc-helix exists in solution or is a crystallization 
packing artifact Notably, this helical region is not observed in the published structures 
of hGHbp (begins at residue 32; deVos et al., Science 255:306, 1992), PRLR (begins at 
residue 2, without any defined secondary structure until the first p-strand, residue 6; 
Somers et al., Nature 372:478, 1994), the INF-yRa (begins at residue 17; Walter et al., 
Nature 376:230, 1995) or the tissue fector (begins at residue 3 without any defined 
secondary structure until the first p-strand, residue 11; Muller et al.. Nature 370:662, 
1994). 

The EBP monomer folds into two domains, Dl and D2, that form an L-shape with the 
long axis of each domain aligned at approximately 90° to each other, the overall 
molecular dimensions are 45 A x 52 A x 62 A. The N-terminal domain (Dl, residues 10- 
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1 14) and C-terminal domain (D2, residues 1 19-220) are connected by a short four residue 
ohdix linker. Both domains are more closely related in overall topology to Fibronectin 
type-in (FBN-m) domains than to Ig domains (Bork et al., J. Mol. Biol., 242:309, 1994). 
The FBN-lll fold is composed of two antiparallel P-pleated sheets, consisting of strands 
A,B E and strands G, F, C and C\ and is found in the two domains of the human growth 
hormone (de Vos et al.. Science 255:306. 1992) and prolactin (Somers et al.. Nature 
372 478, 1994) receptors, the Dl and D2 domains of the a chain of interferon- y receptor 
dFN-YRa) (Walter et al., Nature 376:230, 1995), the D2 domain of CD4 (Wang et al, 
Nature348:411, 1990; Ryu et al., Nature 348:419, 1990), the two domains of tissue factor 
(MuUer et al., Biochemistry 33:10864, 1994; Harlos et al.. Nature 370:662, 1994), the 
third fibronectin-type repeat of tenacin (Leagy et al., Science 258:987. 1992) and the D2 
domain of the chaperone protein PapD (Holmgren et al.. Nature 342:248. 1989). The 
FBN-m topology differs from an Ig constant domain by a shift of strand D from one P- 
sheet (strands A, B, E and D) to the other (strands G, F, C, O. where it is denned as the 
C strand Superposition of equivalent P-sheet core residues of the Dl and D2 domams 
in EBP gives an r.m.s. deviation of 2.3 A for 77 Cec pairs, which is significantly larger 
thanthcc^espor^^ andPRLR(0.8A), and reflects 

- a difference in the subclass of fold between the two EBP domains. 

a- short^h«restdues 10*0), precedes the first P-sandwich that is better 
described as a hybrid of the FBN-m fold with an Ig fold (residues 24-1 14), rather than 
strict FBN-lll topology. In this h-type fold (Wang et al.. Nature 348:41 1. 1990; Ryu et 
aL. Nature 348:419. 1990), the C strand is long and interacts first with strand C and then 
switches to interact with strand E (where C changes Hs designation to strand D) forming 
a four-on-four strand P-sandwich. m^conmr^mP7^r^^^^:bnde^ 

cjg^C^-XP'A)^ 

( between the cysteine pairs that-forrn^ 

compared to 9 and 10 in both GHR and PRLR. The longer connection between strands 
C and E enable the second half of strand C to become strand D. This h-type topology 
is not found in either of the two s-type GHR domains. A potential glycoylation she exists 
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on residue Asn 52 which is located towards the end of the loop region connecting the pB 
and PC strands. Although Asn 52 is not glycosylated in this bacterially expressed protein, 
an external cavity around the Asn 52 side chain could easily accommodate a carbohydrate 
moiety. 

A helical linker (residues 1 15-118) connects Dl to D2 (The <f>, V torsion angles for the 
interdomain helical linker for He'", Asn"\ GIu 117 and Val w are -50° -27°, -76° -21 ° 
99% 26°, and -151 °, 38° respectively.) and has been observed in other members of this 
receptor family, hGHbp, PRLR, IFN-yRa and tissue factor. In EBP, tfce domain 
association is further restricted by a mixed assortment of hydrogen bonding, hydrophobic 
interactions and one salt bridge (between Arg" and Asp 122 ) from 1 1 residues of Dl and 
12 residues of D2 with a total buried surface [The molecular surface areas buried by 
interaction were calculated using the program MS (ConnoDy, J. Appl, Crystallog, 16:439, 
1983) using a 1.7A probe sphere and standard atomic radii (as described in Davies, et al, 
Ann. Rev. Biochem. 59:439, 1990). There may be some discrepancies between values 
reported here and other (deVos et al., Scince 255:306, 1992) published values due to use 
of a different algorithm (Connolly) vs. Lee et aL, J. Mol. Biol., 55: 379, 1971) and probe 
radii. For clarity all values reported here have been calculated in the same way for better 
comparison between the receptors] of 950 A 2 for the two domains. 

D2 (residues 1 19-220) folds into the standard FBN-Ul (s-type) topology with one free 
cysteine and no disulfide bridges, consistent with GHR and PRLR that have three and two 
disulfide bridges, respectively, in Dl but none in D2. After the a-helix linker, D2 begins 
with an irregular coil (residues 1 18-126) that contains Pro 124 which is conserved in the 
structures of hGHbp, PRLR, tissue factor and IFNy-Rcr, and based on sequence 
alignment, in most class- 1 and class-2 cytokine receptors (Bazan, Proc. Natl. Acad. Sd. 
USA 87:6934, 1990). This short coil ends with CHy 124 which has a positive <f> (4>, T * 
52°,40°) consistent with the equivalent Ala 136 and Ala 101 torsion angles in hGHbp (<J>, Y 
c 63°,68°) and PRLR (<fc Y - 58° ; 38°). The Pro" 4 region forms an analogous extended 
bulge conformation adjacent and parallel to a corresponding bulge containing the 



16:33:01 



WO 97/41526 



PCTAJS97/07218 



-15- 



10 



15 



20 



WSXWS motif. The WSAWS sequence forms a modified wide P-bulge (Richardson, 
Adv Prot Chem. 34:167. 1981) and is located in an extended chain region immediately 
preceding the PG strand that would normally connect to the membrane spanning region 
oftheEPOR- 

The quaternary structure of the complex is composed of two peptides and two receptors 
that form a T-shapes assembly. A noncovalent peptide dimer interacts with two receptor 
molecules to generate an almost perfect 2-fold symmetrical arrangement. After 
superposition of D2 of the two EBP molecules in the dimer. the centers of mass of the 
two Dl domains are only 0.8 A apart, sufficient to perturb perfect two-fold symmetry. 
Separate superposition of the corresponding Dl and D2 of each receptor in the dimer 
results in r.m.s. deviations of 0.53 A (105TM C« pairs) and 0.47 A (93 D2 Ca pairs). 

The cyclic EMP1 contains a single disulfide bridge between Cys" and Cys™. which links 
two short P-strands (residues 4-7 and 13-16) that are connected by a slightly darted 
type 1 P-tum [Pro- 0 (M) and Leu" 1 <M) of the P-turn have 4».<P * -62% -38« and -99 . 
-60<> respectively. The carbonyl oxygen of Leu'" has a hydrogen bond to EBP distorting 
the T value from its normal 0«*30- (H-2) in a standard type I P-turn.] consistmg of 
residues Gly w -Pro"°-Leu" , -Thr m . Each peptide has a very close association with its 
other peptide partner and buries 320 A of its 1220 A 1 molecular surfcee in this interaction 
(Connelly, J. Appl. Crystallog. 16:439, 1983; Davies et aL, Ann. Rev. Biocherr, 59:439, 
1990- Richards, J. MoL Biol. 55:379, 1971). Four hydrogen bonds between the 
^chains of the two peptides results in formation of a four-stranded anti-parallel P- 
pleated sheet (Table 2). Two symmetric hyrdopbobic cores are assembled by peptide 
dimerization and are comprised of the disulfide bridges and the side chains of Tyr" J»he 
and Trp ro The construction of each hydrophobic core resembles a box which places the 
aromatic rings of Phe".Trp™ and Tyr~ (from the other peptide) and the disulfide bridge 
(Cys M -Cys" J ) at the corners. The two glycine residues at either end of the pept,de are 
not structured. 
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Thc peptide dimer is embedded in a deep crevice between two EBP receptor molecules 
A portion of each peptide monomer interacts with both receptor molecules. The binding 
sites of each EBP are practically identical due to the 2-fold symmetric interactions 
imposed on binding the peptide dimer. The four major contact areas on EBP come from 
segments on four loop regions (LI. L3, L5. L6) that connect strands A to B (LI residues 
33-34) and F to G (L6 residues 90-94) in Dl and strands B to C (L5 residues 148-153) 
and F to G (L6 residues 203-205) in D2. The total buried molecular surfaces in the 
peptide-EBP assembly are 840 A 1 and 880 A 2 for the two peptides and EBP*s. respec- 
tively. The peptide-EBP interaction can be separated into distinct hydrophobic (67%) and 
polar (33%) areas. A hydrophobic core is formed between the peptide and receptor and 
comprises Phe",Met 150 and Phe 1 " from one EBP molecule and the peptide hydrophobic 
box consisting of Phe M and Trp WJ from one peptide and Tyr" and Cys" 5 from the other 
peptide. The polar interactions are located mainly at the bottom of the binding crevice 
and are mainly with loop L5 in D2. Five of the six hydrogen bonds are between the 
mainchain of the P-turn residues Gly w ,Pro ,,,0 and Led" from one peptide with the 
mainchain and sidechain hydroxyl of conserved Tyr™, which crosses over its other 
peptide partner, to interact with loop L3 (Table 2). The EBP-EBP interaction makes a 
surprisingly minor contribution to the overall stability of the complex where the inter- 
receptor buried molecular surface is only 75 A*, contributed by Leu ,7J and Arg"' from 
each receptor molecule. 



EMP1 is one of a femih/ of sequences that contain several conserved residues, besides the 
c y steines <^***»EftMM«MM*^ ^ The most structurally signifies! 

of these consensus residues appear to be Tyr* 4 and Trp™ which along with the disukffl 
bridge have a major contribution to the hydrophobic core of the peptide-peptffi! 
interaction. Moreover, these two aromatic residues play a pivotal role in peptide-recepRr 
interaction and in receptor dimerization. 



Dimerization of EBP in Solution 
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To explore the interaction of EMP1 with EBP in solution we employed a 
bifunctionalsulphydryl reactive crosslink* DPDPB, [l,4-di-<2^pyridyldithio 
propionamido) butane], in an attempt to stabilize a peptide-dependent dimeric structure. 
The choice of crosslinker was based on previous experiments with amine-reactive 
crosslinkers that were found to inactivate EBP. EBP contains a single free sulphydryl 
(Cys m ) in D2 which is potentially reactive to crosslinking reagents (The DPDPB 
crosslinker itself does not inactivate the EPO binding potential of EBP nor the 
proliferative properties of EMP1). A dimeric EBP product is formed by co-incubation 
of EMP1, DPDPB and EBP. The amount of dimeric product increases with peptide 
concentration and no significant dimer product is observed in the absence of peptide. 
DPDPB-crosslinked products formed through disulfide-exchange reactions should be 
readily reversible by reduction as is seen for the covalently-linked EMPl-mediated dimer. 
Furthermore, we have constructed a covalently-linked dimeric form of EMP1 that 
demonstrates increased biological potency (Johnson et al. in preparation). The Cys 1 " 
residues in D2 of the EBP dimer are 20.7 A apart (S Y -S Y distance) which approximates 
the 16 A length (and approximately 2 A in bond length at each end) of the DPDPB cross- 
linker. Thus EMP1 mediates formation of a soluble EBP dimer complex in solution 
consistent with the crystal structure. 

TheWSXWS motif 

The WSAWS sequence (residues 209-213) corresponding to the WSXWS box occurs in 
a P-bulge (Richardson, Adv. ProL Chem. 34:167, 1981; Chan et al., Protein Science, 
2: 1 574, 1993) immediately preceding P-strand G in D2. Residues in this motif do not 
interact with Bgand, have no role in receptor-receptor interactions and are located on the 
opposite side of the receptor-receptor and receptor-ligand interfece. The WSAWS box 
represents only a segment of a complex array of interactions that involves several other 
conserved side chains from the four-stranded P-sheet in D2. The indole ring systems of 
Trp™ and Trp"* point toward an external concave surface of the p-sheet and are only 
partially solvent exposed, whereas the Ala 111 side chain points directly out into solution. 
The amides and hydroxyls of both Ser™ and Ser™ form hydrogen bonds with the main 
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chain of residues 198 and 196 of adjacent strand F in a pseudo p-sheet type interaction 
that resembles a modified wide p-bulge (Richardson, Adv. Prot. Chem. 34:167, 1981; 
Chan et al., Protein Science, 2:1574, 1993) where the sidechain hydroxyl rather than the 
carbonyl oxygen makes the p-sheet interaction. The P-budge architecture places the two 
Trp residues, which are spread four residues apart, on the same side of the p-sheet and 
not on opposite sides as in normal p-sheet or extended chain structures. The guanidinum 
group of Arg 157 from Strand F, the central residue (Richardson, Adv. Prot. Chem. 34: 1 67, 
1981; Chan et al., Protein Science, 2:1574, 1993) in the bulge, is positioned exactly 
between the two Trp indole rings to form an extended 7c-cation system (Kumpf et al., 
Science 261:1708, 1993. The center of the pyrrole ring of Trp 209 , the Ne of the Arg 197 and 
the center of the benzene ring of Trp 2u are positioned on a straight line with the three 
planes of the conjugated systems stacked parallel to each other at approximately 4 A 
spacing. In addition, the aliphatic portion of the Arg 199 side chain has hydrophobic 
interactions with the indole ring of Trp 209 , completing the alternating stacking of two 
aromatic and two positively-charged amino add residues. The side chain of Glu 157 forms 
a hydrogen bond with Arg 197 presumably to help orient the guanidinium group and add 
some specificity and stabilization to the system. 

It appears then that the linear WSXWS motif identified from sequence alignments of 
cytokine receptors represents only a component of a more complex conformational unit 
that contributes a significant structural feature to D2. Aromatic residues have previously 
been suggested to have a stabilizing effect and play a role as a folding nuclei in structures 
of anti parallel p-sandwiches (Finkdstein et al., Protein Eng. 6:367, 1993). The amino- 
aromatic parallel stacking between the guanidinium group of arginine and the aromatic 
rings is a common feature in protein structures (Burley et al., Adv. Prot. Chart, 39: 125, 
1988; Flocco et aL, J. MoL BioL, 235:709, 1994), but a parallel triple stacking of n-cation 
systems is rare (Kim et al., Biochemistry 32:8465, 1993) although observed in other 
class- 1 cytokine receptors, hGHbp and PRLR. 
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The structural equivalents of the WSXWS motif in hGHbp (YGEFS) and PRLR 
(WSAWS) are involved in an even more intricate and complex array of rc-cation 
interactions. The it-cation system is extended in hGHbp and PRLR to include an 
additional aromatic residue <Jn> tu for hGHbp and Trp ISS for PRLR) from the loop region 
that links PC and PC in D2 and a positively-charged residue (Arg 211 for hGHbp and 
Arg 147 for PRLR) that stacks between the Trp and the second aromatic residue. The 
additional Arg residue is contributed either from the pF strand as in hGHbp (Arg 211 ) or 
from PC as in PRLR. (Arg 147 ); the glutamine residue that hydrogen bonds and orients the 
arginine also switches strands. Sequence alignments suggest that this Arg-Gln switch 
could be common to other members of the class- 1 cytokine receptor family. The 
extended n-cation system in hGHbp and PRLR consists of five positively charged and 
three aromatic residues stacked in an alternating order which comprises of Lys 215 , Tyr 222 , 
Arg 1 ", Phe 221 . Arg 2, \ Trp 1 **, Lys' 7 ' for hGHbp and Lys" 5 , Trp ,, \ Arg ,n , Trp"\ Arg 147 , 
Trp ,ss , Lys 14 ' for PRLR- The first aromatic- Arg-aromatic trio are approximately 4 A apart, 
as in EBP, but the second system is stacked closer together at approximately 3.6 A 
spacings consistent with it-* interaction (Burley et al., Adv. Prot. Chem., 39:125, 1988; 
Flocco et al., J. Mol. Biol., 235:709, 1994). The outer lysines also use the aliphatic 
portions of their side chains to form hydrophobic interactions whh the aromatic rings. 
Based on sequence alignments with other members of the class- 1 cytokine receptor 
superfamHy, such structurally extended u-cation systems could exist in human 
thrombopietin, EL-6 and ciliary neurotrophic factor receptors, and in human IL-4 receptor 
based on structural modeling (Gustchina et al.. Proteins 21:140, 1995). Although EFN- 
YRa and tissue factor do not have a WSXWS motif, the corresponding sequences TTEKS 
(residues 213-217) for IFN-YR« (Walter et al.. Nature 376.230, 1995) and KSTDS 
(residues 201-205) for tissue factor (MuUer et al., Biochemistry 33: 10864, 1994; Harlos 
et al., Nature 370:662, 1994), maintain a very similar p-bulge. The consensus sequence 
among these five x-ray structures indicates that a serine or threonine in positions 2 and 
5 maintain a common set of hydrogen bonds between their side chain hydroxyls and the 
mainchain of the neighboring strand. Only in hGHbp is there no hydroxyl-containing 
residue in position 2, but Ser 226 still maintains the equivalent interaction. A Ser 22 ' to Ala 
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mutation abrogates hGHR binding to hGH, and its expression on the cell surface is 
drastically reduced (Baumgartner et al. t J. Biol. 'Chan., 269: 29094, 1994). In GM- 
CSFRa and IL-2R0, point mutations of the serine residues cause a substantial decrease 
in cell surface expression but little or no effect on ligand binding (Ronco et al, J. Biol 
Chem. 269:277, 1994; Miyazaki et al., EMBO Journal 10:3191, 1991). 

Conservation of the WSXWS motif in EPOR or its equivalent in other members of the 
class 1 cytokine receptors has been proposed to be essential for biological activity and 
was thus assumed to be part of the receptor binding site CYoshimura et al., J. Biol. Chem. 
267:1 1619, 1992; Quelle, MoL Cell. Biol. 12:4553 1992). For EPOR, a systematic study 
of 100 mutations of the WSAWS sequence demonstrates that most of the mutations of the 
two tryptophan and serine resulted in molecules that did not reach the cell surface but 
were retained in the endoplasmatic reticulum (Hilton et al., Proc. Natl. Acad. Sci. USA 
92:190, 1995; Hilton et al., J. Biol. Chem. 271:4699, 1996). Furthermore, an Ala 211 to 
Glu mutation in the WSAWS sequence resulted in better transportation from the ER to 
the Golgi and a 3-5 fold increase of the number of EPOR molecules on the ceil surface 
compared to the wild-type (Hilton et al., Proc. Natl. Acad. Sci. USA 92: 190, 1995; Hilton 
et al., J. BioL Chem. 271:4699, 1996). These results support our conclusion that the 
WSXWS sequence plays an important role in the structure and folding of D2 in EPOR 
and other related receptors. 

Comparison with other cytokine-receptor complex structures 
The overall quaternary structure of the peptide-EBP complex substantially from the 
equivalent arrangement in the hGH-hGHR complex. The non-symmetric nature of the 
single four-helix-bundle structure of the growth hormone ligand results in an asymmetric 
homo-dimerization of the receptor that corresponds to a 159° rotation between receptors 
compared to the almost perfect 2-fo!d (1 80°) rotation for the EBP-peptide complex. The 
tertiary arrangement of domains within EBP and hGHbp is also somewhat different. 
When the equivalent EBP and hGHbp D2 domains are superimposed on each other, their 
corresponding DI domains differ by a 12° rotation and a 4.3A translati n. 
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The mechanism of hGH binding to its receptor has been well studied (Wells, Curr. Opin. 
Cell Biol. 6:163. 1994; Clackson et al., Science 267:383, 1995) and is sequential. Initial 
high affinity (nM) binding of the hormone with one receptor results is a buried surface 
of 1 130 A 1 on the receptor. The second hGHb P 2 has a substantially smaller interface 
(deVos et al.. Science 255:306, 1992) with the second binding site on hGH and interacts 
only with the preformed 1:1 complex to generate buried surface areas of 740 A J with hGH 
and 440 A 1 with the first hGHbpl (deVos et al.. Science 255:306, 1992; (Wells, Curr. 
Opin. Cell Biol. 6:163. 1994; Clackson et al.. Science 267:383. 1995). The binding 
determinants of each hGHbp are comprised of the six recognition loops (L1-L6). three 
of which (L1-L3) come from one end of the P-sandwich structure in Dl. one from the 
interdomain linker and two from D2. 

Although these two receptor complexes. EBP-EMP1 and hGH-hGHbp, have different 
dimeric arrangements, which probably in this case represent differences in the size and 
shape of the natural versus synthetic ligand. both receptors share equivalent ligand 
recognitionloops.Ll.L3.L5andL6fortheEBPandLltoL6forthehGHbp. Anon- 
active PRLR, complexed with only one molecule of hGH, also uses the same contact 
loops (LI to L6) (Somers et aL, Nature 372:478, 1994)7 Based on similarity of the ligand 
recognition sites in hGHbp and PRLR, one would expect that the binding site of EBP, 
when its natural EPO ligand is bound, would extend to include two additional loops, L2 
and L4. mat comprise residues 59-63 (L2) between strands C to C. and residues 110-118 
(L4) from the carboxyl end of PG in Dl and the interdomain linker. These six loops in 
EBP, hGHbp and PRLR area in structurally equivalent positions but vary in size, amino 
acid composition and conformation although the interacting portions of each loop (side 
or tip) remain similar, LI, L2, L3. L5 interact mainly with their tips and L6 with its side. 
In EBP, the L5 loop is three residues shorter than in hGHbp and PRLR, where the L6 
loop is three and four residues longer than in hGHbp and PRLR, respectively. The L2 
loop also varies (6 to 10 residues) among the three receptors but in EBP does not 
participate in peptide binding, and in hGHbp is partially disordered, although it does 
contact th hormone. In one respect, this situation is similar to the complementarity- 
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determining regions (CDR's) in antibodies, where changes in length and sequence of the 
six binding loops impose specificity for different antigens, whereas the framework itself 
remains constant (Wilson et al. t Ciba Foundation Symposium. Wiley, Chichester, 1991, 
Vol. 159, p. 13). 

It has been shown for the hGH-hGHbp complex that only a subset of 9 out of 33 
interacting residues that make up the structural epitope of the receptor constitute a 
functional epitope or hot spot (Wells, Curr. Opin. Cell Biol. 6: 163, 1994; Clackson et al., 
Science 267:383, 1995) where high affinity binding interaction takes place. This reduced 
epitope is substantially smaller than the structural epitope and is comprised from residues 
(Arg* Glu", lie 103 , Trp I0 \ lie 105 , Pro 106 . Asp 163 , and Trp 1 ") which are located in contact 
loops LI, L3 and L5 with the most significant contribution (>4.5 kcal/mol) coming from 
two aromatic residues (Trp 104 and Tq* 69 ) in L3 and L5 (Wells, Curr. Opin. Cell Biol. 
6:163, 1994; Clackson et al. # Science 267:383, 1995; Wdls, Proc. Natl Acad. Sci. USA 
93:1, 1996). In EBP, Phe 93 is equivalent to Trp 104 in hGHbp, as suggested previously 
(Wells, Cunr. Opin. Cell Biol. 6:163, 1994; Clackson et al., Science 267:383, 1995; Wells, 
Proc. NaU Acad. ScL USA 93:1, 1996; Jolliflfe et al., Nephrol Dial. Trans. 10:suppl. 2, 
28, 1995), but there is no homologous residue to Trp 1 ** in the shorter L5 loop. In the 
EBP-EMP1 complex, the Phe M peptide aromatic side chain occupies the equivalent 
position of the Trp 1 * 9 side chain in hGHbp. One can assume that when EPO binds to its 
receptor, the hormone may provide an aromatic residue to the hydrophobic core of the 
binding interface and/or the L6 loop in EBP may play a more significant role in the 
hormone binding than in hGHbp, since it is 3 residues longer and contains the aromatic 
Phe 205 . 

In these three class- 1 receptor structures, some loops are disordered which are in D2 for 
EBP for EBP (residues 164-166 in EBP1 and 133-135 in EBP2) and in Dl for both 
hGHbp (residues 55-58, 73-78 for hGHbp 1 and 54-60, 73-75 for hGHbp2) and PRLR 
(residues 31-33, 84-86). Otherwise, these three class- 1 cytokine receptors do not differ 
greatly in their over all tertiary structures; Dl and D2 have broadly similar general 



16:33:01 



• 



PCT/US97/07218 

WO 97/415Z6 



-23- 



10 



15 



20 



arrangement in all three receptors such that the angle between the long axes of the two 
domains is approximately 90 degrees. I tis this arrangement of domains that allow these 
particular L1-L6 loops tobe available for the recognition and binding of ligands. In a 2:2 
complex between IFN-Y and its class-2 receptor IFN- Y R«, ™ and ™ are related by a 
125 degree angle, which elongates the receptor and restricts the binding determinants that 
can be used for interaction with hormone; the LI loop now becomes buried in the D 1 -D2 
interface although the other five loops (L2-L6) are still available for ligand interacts 
This elongated interdomain arrangement is also observed in tissue factor (MuUer et al.. 
Biochemistry 33:10864. 1994; Harlos et al.. Nature 370:662. 1994) which has a distant 
relationship to the cytokine receptor superfamily. 

A mutational analysis of the EBP molecule indicates that the most crucial amino acid 
residue for binding EPO is Phe» in the L3 loop (Jolliffe et aL. Nephrol. Dial. Trans. 
lOsuppl 2.28. 1995). The Phe93Ala mutant shows an increase int he IC* compared to 
the wild-type by a factor of approximately 1000. whereas other mutants (Ser91Ala. 
Scr92Ala, Val94Ala, Metl50Ala and Hisl53Ala) show small relative increases m teh 
ICo of only 2.5-12.5 fold). The side chain of Phe« buries 66 X of molecular surface, 
which is the highest among interacting side chains. In hGHbp. the corresponding 
TrplWAla mutation results in an increase in the K, by a factor of more than 2.500 
compared to the wild-type indicating the equivalent importance of this residue in hGH 
binding and its key contribution to the hydrophobic core of the functional epitope (Wells, 
Curr. Opin. Cell Biol. 6:163. 1994; Clackson et al.. Science 267:383. 1995; Bass et al. 
Proc. Natl. Acad. ScL USA 88:4498. 1991). 

The role of dimerization on signal transduction 

In the EBP-EMP1 complex structure, we surprisingly observe that a peptide, unrelated 
i„ sequence and probably m structure, to the natural ligand. can induce a biology 
active dimerization of EPO receptor that promotes signal transduction and cell 
proliferation. Comparison of three dass-1 cytokine receptor complexes, whose structures 
hav been determined so far. suggests that when the natural EPO hormone, wluch .s 
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proposed to have a structure of a four-helix bundle (Boissel et al., J. Biol. Chem. 
268:15983, 1993), induces receptor dimerization, it is more likely to resemble the hGH- 
hGHbp assemblage. This would suggest that more than one mode of productive 
extracellular dimerization is permissive for intracellular dimerization of the cytoplasmic 
domains with two JAK2 molecules in order to initialize the cascade of events that 
produces the biologically relevant signal (Ihle et al., Seminars in Immunology 5:375, 
1993; Klingmuller et al., Cell 80:729. 1995). The peptide-EBP structure would then' 
represent only one possible dimeric arrangement that promotes signal transduction. 

Mutant EPOR molecules, containing a single Arg to Cys mutation (Arg" 0 in human and 
Arg 1 * 9 in murine), have been shown to form biologically active dimers in the absence of 
EPO (Yoshimura et al., J. Biol. Chem. 267:11619, 1992); Watowich et al.. Proc. Natl. 
Acad. Sci. USA 89:2140, 1992; Watowich et al., Mol. Cell. Biol. 14:3535, 1994). 
suggesting that extracellular recptor homo-dimerization may be sufficient in itself for 
signal transduction. It has been shown in another system (Spencer et al.. Science 
262:1019, 1993) that activation of a specific set of trasncription factors can be induced 
by the chemical crosslinking of cytoplasmic domains of modified cell membrane 
receptors that do not contain the extracellular and transmembrane domains. These 
receptors are not related to the cytokine receptor superfamily but illustrate that 
oEgomerication plays a key role in activation of the receptor, and that the main functional 
role of the extracellular, ligand-binding domain is to allow Cm the presence of ligand) 
dimerization or oligomerization and induce similar association of the cytoplasmic 
domains. 

Mutageneses experiments originally suggested a role for the WSXWS motif in this cell 
signalEngr»ocess (Yoshimura etal., J. Biol. Chem. 267:11619, 1992; Quelle et aL MoL 
Cell. Biol. 12:4553, 1992; Chiba et al., Biochera. Biophys. Res. Comm. 184:485, 1992) 
possibly by promoting receptor homo-dimerization. However, truncation mutants of 
EPOR (Miura et aL, Arch. Biochem. Biophys. 306:200, 1993) do not confirm this role for 
the WSXWS motif. The EBP-EMP1 complex structure shows that the WSXWS motif 
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of the EPOR, as for the hGH-hGHbp complex (deVos et al., Science 255:306, 1992) is 
located on the opposite face of the molecule from the receptor dimerization. In the 
absence of unliganded structures for the extracellular domains of EPOR, hGHR and 
PRLR, it is not possible to determine whether any conformation change occurs on ligand 
binding that would involve the WSXWS box. Apart from being a striking structural 
feature in D2, and its obvious proximity to the membrane spanning domain, one cannot 
rule out possible interactions of this region with some other cell surface molecules that 
are involved somehow in the signal transduction process. 

Towards design of small molecule mimetics 

*the*«itura^ indu« the appropriate biological response.' 

The peptide can be assumed to form a substantially smaller contact interface than the 
natural hormone with the receptor. The peptide binding site in EBP forms an almost flat 
surface which is mainly hydrophobic in nature, without any cavities or charged residues 
that are normally essential for the specific targeting of small molecule ligands to a 
receptor binding site. The hGHbp study (Wells et al.. Science 267:383, 1995; Wells. 
Proc. NatL Acad. Sci. USA 93:1, 1996) shows that only a small part ; of the observed 
structural binding site, the so-called functional epitope (SUETS), contributes most of the 
binding energy and strongly implied that a -nnnimized- hormone designed to interact 
with this site could form sufficient interactions to activate the receptor. Furthermore, the 
Umited site of interaction of the small agonist peptide with the EBP corresponds almost 
exactly to the smaller functional epitope derived from alanine scanning of hGH and 
hGHbp. fb**?****^^ we lu.ve arrived at the similar conclusion that a 
smalhnumber of key interactions can contribute to a functional epitope on a receptor. 
Understanding of this simplified interaction surface can be now combined with further 
mutational studies to assist in identifying the most crucial residues in the functional 
epitope, and consequently provide a more practical target for drug design. 

Data Collection, MER and Refinement Statistics 



16:33:01 



WO 97/41526 



PCT/US97/07218 



-26- 



Thc crystallographic data is summarized in Table I. Native crystallographic data were 
collected on a Siemens multiwire area detector mounted on an Elliott GX-18 generator 
operating at 40kV and 55mA, with a crystal-to-detector distance of I20mm. Two 
derivative data sets were collected on a MAR image plate mounted on a Siemens 
generator operating at 50kV and 80mA, with ciystal-to-image plate distance of 1 50mm. 
Data were integrated, scaled and reduced using the programs XENGEN (Howard et al. 
J. App. Cryst. 20:383, 1987) for the native data and DENZO/SCALEPACK (Otwinowski 
et al, SERC Darsbury Laboratory, Warrington, 1993) for the derivative data. Initial 
multiple isomorphous replacement anomalous scattering (MIRAS) phases were calculated 
to 3.1 A using the program package PHASES (Furey, American Crystallographic 
Association Fortieth Anniversary Meeting, New Orleans, LA, 1990) with a mean figure 
of merit of 0.64 (25.0-3.1 A). Phases were refined in PHASES using the solvent 
flattening protocol to a mean figure of merit of 0.92 (25.0-3.1 A). The quality of the map 
was generally good and most of the complex structure (94%) could be fitted using the 
graphics program O (Jone et ai., Acta Crystallogr A47:l 10, 1991). The register of the 
amino acid residues was verified from the positions of the two disulfide bridges in Dl, 
and the positions of the two Hg"s from the mercury acetate derivative that were correctly 
assumed to bind to the free Cys 1 " residue, the peptide interpretation was verified from 
another data set from a complex between EBP and an iodinated peptide (Tyr p4 was 
substituted for/Modo-PheX which diffracted to 3. 3 A resolution, that in difference Fourier 
(Fbdo-F^a,^^ gave a clear indication of the location of the iodine atoms. The structure 
was refined using the slow-cooling protocol in X-PLOR 3_l(Bmnger et aL. Acta 
Crystallogr A46:585, 1990; Brunger, X-PLOR, Version 3.1: A System for X-ray and 
NMR, Yale Univ. Press, New Haven, CT, 1992) and rebuilt using Fo-Fc, 3Fo-2Fc and 
SIGMAA(Read, Acta Crystallogr. A42:140, 1986) weighted electron density maps. After 
every two cycles of refinement, a set of simulated annealing omit maps (7-10%) to reduce 
model bias was calculated and the entire structure rebuilt.. After several cycles of 
refinement, individual temperature factors were calculated and after 10 cycles of 
refinement and model building, the R-value was 0.21 for 8.0-2.8 A data with F>lo 
(13,984 reflections). The average thermal parameters for receptor I, receptor H and the 
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peptides are 10.5A 2 , 12.3A and 10.7* respectively. Only one non-glycine residue 
[Asn 1 * 4 in EBP2], located in a loop region in Dl, is in a disallowed region in the 
Ramachandran plot. No solvent molecules were included in the model due to the 
moderate resolution (2.8 A) of the structure determination. 

Binding Contacts 

Binding contacts are summarized, in part, in Table 2: Hydrogen bond interactions in the 
binding she of the EBP-EMP1 complex. Due to the symmetrical nature of the complex, 
peptide-1 and peptide-2 have equivalent interactions with the two EBP molecules. The 
hydrogen bond interactions were analyzed using HBPLUS (McDonald et al.. J. Mol. Biol. 
238:777, 1994), based upon both distance (3.9 A cutoff) and geometrical considerations. 



A number of embodiments of the present invention have been described. Nevertheless, 
h will be understood that various modifications may be made without departing from the 
spirit and scope of the invention. Accordingly, it is to be understood that the invention 
is not to be limited by the specific illustrated embodiment, but only by the scope of the 
15 appended claims. 
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4 PHE 308 CE1 1 3.95 

4 PHE 308 CD1 1 4.11 

4 TRP 313 CH2 1 3*74 

4 TRP 313 CH2 1 3. 98 

4 TRP 313 CZ2 1 4 0 8 

4 GLY 309 O 1 3*59 

4 GLY 309 O 1 3*49 

4 GLY 309 O 1 3 * 66 

4 PRO 310 O l 3 35 

4 PRO 310 C 1 3*62 

4 LEU 311 O 1 3*41 

4 GLY 309 O 1 3*78 

4 LEU 311 C 1 3*87 

4 PHE 308 CD2 1 3*50 

4 PHE 308 CB 1 3 "70 

4 PHE 308 CG 1 3*79 

4 PHE 308 CD2 1 3*52 

4 THR 312 C 1 3*55 
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4 PHE 308 CD2 1 3*45 
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4 PRO 310 0 1 3*82 
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A LEU 311 O 1 3 54 

4 LEU 311 O 1 3*57 
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4 THR 312 CB 1 3. 4S 

4 THR 312 CA 1 3.76 
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4 PHE 308 CE2 1 3*40 

4 PHE 308 CZ 1 3.53 

4 TYR 304 OH 1 3.44 

4 TYR 304 OH 1 3.88 

4 PRO 317 CB 1 3 95 

4 TYR 304 OH 1 . 3*44 

4 PRO 317 CB 1 3.6I 

4 TYR 304 CZ 1 3.83 

4 TYR 304 CE2 1 3.34 

4 TYR 304 OH 1 3 62 

4 TYR 304 CE2 1 3 66 

4 TYR 304 CZ 1 3-e8 

4 TYR 304 OH l 3. 80 

4 TYR 304 OH 1 3.73- 

4 TYR 304 CE2 1 4. 00 

4 TYR 304 CE2 1 3 53 

4 PRO 317 CD 1 3*59 

4 CYS 315 O 1 3*74 

4 CYS 315 CB l 3^55 

4 TYR 304 CD 2 l 3^72 
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CLAIMS 

What is claimed is: 

1. A computer-assisted method for identifying potential mirnetics of erythropoietin, 
using a programmed computer comprising a processor, a data storage system, an 
input device, and an output device, comprising the steps of: 
(a) inputting into the programmed computer through said input device data 
5 comprising the three-dimensional coordinates of a subset of the atoms in 

the peptide GGTYSCHFGPLTWVCKPQGG when said peptide is co- 
crystallized with a portion of the erythropoietin receptor comprising 
amino adds 1 to 225 of said receptor, thereby generating a criteria data 
set; 

10 (b) comparing, using said processor, said criteria data set to a computer 

database of chemical structures stored in said computer data storage 
system; 

(c) selecting from said database; using computer methods, chemical structures 
having a portion that is structurally similar to said criteria data set; 
1 5 (d) outputting to said output device the selected chemical structures having 

a portion similar to said criteria data set. 
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2. A computer-assisted method for identifying potential mimetics of erythropoietin, 
using a programmed computer comprising a processor, a data storage system, an 
input device, and an output device, comprising the steps of: 
(a) inputting into the programmed computer through said input device data 
5 comprising the three-dimensional coordinates of a subset of the atoms in 

the peptide GGTYSCHFGPLTWVCKPQGG when said peptide is co- 
crystallized with a portion of the erythropoietin receptor comprising 
amino adds 1 to 225 of said receptor, thereby generating a criteria data 
set; 

10 0>) constructing, using computer methods, a model of a chemical structure 

having a portion that is structurally similar to said criteria data set; 
(c) outputting to said output device the constructed model. 

3. A compound having a chemical structure selected using the method of claim 1, 
1 5 said compound being an EPO mimetic. 

4. The compound of claim 3 wherein said compound is not a peptide. 
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6. 



The compound of claim 3 wherein said compound is a peptide. 

The compound of claims 5 wherein said peptide has 15 of fewer amino acids 
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